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We investigate thermal and magnetic properties of the double spin chain com- 
pound KCuCls via an exactly solved ladder model with strong rung interaction. 
Results from the analysis of the thermodynamic Bethe Ansatz equations suggests 
the critical field values H^i = 22.74 T and Hc2 = 51.34 T, in good agreement with 
the experimental observations. The temperature dependent magnetic properties 
are directly evaluated from the exact free energy. Good overall agreement is seen 
between the theoretical and experimental susceptibility curves. Our results sug- 
gest that this compound lies in the strong dimerized phase with an energy gap 
A 35 K at zero temperature. 



1. Introduction 

It is believed that the compounds KCuCls, TlCuCls and NH4CUCI3 exhibit 
a double spin chain structure, ^^^^ along the lines of Fig. 1. In the double 
chain structure, coupling constants J_l {J\\) denote the interchain (intra- 
chain) spin exchange interactions, with Jd a diagonal interaction. However, 
there appears to be no uniform agreement on the values of these coupling 
constants for the double chain compounds. In particular, the coupling 
constants for the compound KCuCls arc uncertain. Several theoretical 
models have been proposed to describe this material, including a double 
chain model with strong antiferromagnetic dimerization,^ a ladder model 
with additional diagonal interactions^'^ and a three-dimensional coupled 
spin-dimer system.^"* None of these models provide an overall fit for all 
thermal and magnetic properties, see, e.g., the review by Dagotto.^^ Mea- 
surements of the high field magnetization^'^ and the susceptibility^ indicate 
that KCUCI3 exhibits a singlet ground state with an energy gap A w 31 K 
at T = 1.7 K. Nevertheless, it has been difficult to fix all of the coupling 
parameters of the model by fitting to only one physical property at a time. 
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Figure 1. Schematic picture of the structure of double chain compounds such as 
KCuCls. Here Jx {J\\) is the interchain (intrachain) interaction. is the spin ex- 
change interaction in the diagonal direction. 

At very low temperatures, T < 5 K, the compound KCuCls exhibits three- 
dimensional magnetic ordering due to complex structural magnetic inter- 
action paths. 

In this communication we investigate the critical fields, magnetization 
and susceptibility of the compound KCuCls via an integrable ladder model. 
The results are used to examine the values of the coupling constants for the 
double chain structure. The results for the ladder model with strong rung 
coupling are seen to be in good agreement with the experimental results 
for the energy gap, critical fields, susceptibility and magnetization. 

2. The integrable ladder model 

It has been shown that integrable (exactly solved) ladder models can be 
used to describe real ladder compounds with strong rung interaction. 
These integrable ladder models enjoy the nice property that thermal and 
magnetic quantities can be obtained exactly via well developed nicith- 
ods from integrable systems, such as the Thermodynamic Bcthc Ansatz 
(TBA),!*^ the Quantum Transfer Matrix (QTM),i^ T-systems^*^ and the 
High Temperature Expansion (HTE) of Non Linear Integral Equations 
(NLIE).i9-22 

The simplest integrable two- leg spin- ^ ladder model is constructed from 
the integrable su(4) spin chain with singlet rung interaction. The Hamil- 
tonian is given by^^ 

L L 

H = J||Hieg + J±Y1 ■ - f^B9Hj2iSj + (1) 

j=l 3=1 

where 

L 

Hieg = J2 {^j ■ ^0+1 + Tj ■ fj+, + 4(5,- • Sj+r){fj . f,+i)) . (2) 

3 = 1 



February 2, 2008 18:47 Proceedings Trim Size: 9in x 6in 



NANKAI 



3 



Here L is the number of rungs with Sj = {S^, S], S?) and fj = (T?, Tf, T?) 
spin-^ operators acting on site j. The Bohr magneton is and g is the 
Lande factor. Periodic boundary conditions, Sl+i = Si, T^+i = Ti, are 
imposed. 

In contrast to the standard Heisenberg ladder model the integrable lad- 
der model features an additional biquadratic spin interaction term in the 
definition (2) of Higg. This term causes a shift in the critical value of the 
rung coupling J± at which the energy gap closes, and it also causes a rescal- 
ing of the parameter J|| for strong rung coupling. In the strong coupling 
limit J± ^ J|| the rung interaction dominates the ground state and low- 
lying excitations. The integrable model then lies in the same phase as the 
standard Heisenberg ladder, motivating its analysis. 

The ground state properties at zero temperature may be obtained from 
the TEA equations. ^'^'■^''"^^ Details of the derivation can be found in Ref. 15. 
In the strong coupling limit the integrable spin-| ladder model exhibits 
three quantum phases: a gapped phase in the regime H < Hci, a fully 
polarized phase for H > Hc2 and a Luttinger liquid magnetic phase in 
the regime Hci < H < Hc2- The exact values for the critical fields are^^ 
Hci = Jj_ - 4J|| and Hc2 = J±+ 4J||. 

On the other hand, the temperature dependent free energy has been 
calculated via the exact HTE of the NLIE.^**'^^ The free energy of the 
integrable spin ladder (1) is given in the form^*'^^ 



(3) 



where Q^^^ and the first few coefficients c^^\ are given explicitly in Refs. 14, 
15. These terms are functions of the rung coupling J±, nsgH and the 
temperature. Most importantly, the exact expression (3) for the free en- 
ergy can be used to examine physical properties such as the magnetization, 
susceptibility and magnetic specific heat via the standard thermodynamic 
relations 



M = - 



df{T,H) 



dH 



x = - 



d^f{T,H) 



C=-T 



H 



3. The compound KCuCla 



In this section we examine the low temperature properties of the compound 
KCuCls- Experimental measurements of the high field magnetization^'^ 
show that magnetic anisotropies are negligible, because the critical fields 
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are almost the same for the external field in different directions. However, 
the susceptibility curves for the external magnetic field along the differ- 
ent directions are influenced by different ^-factors.^ In this way magnetic 
anisotropics may lead to different critical fields for external magnetic fields 
along different directions. This can be easily seen from the TBA analysis. 
For instance, if the rung interaction along the z-axis is increased, i.e., by 
adding an extra term = J2j=i '^j'^j the rung interaction, the critical 
fields for the magnetic field along the ^-direction are given by 



The experimental results^'^'^ indicate that is negligible. Analysis of 
such anisotropic behaviour can be found in Rcf. 26. We therefore take the 
high field magnetization curves for the external field along the perpendicu- 
lar and parallel directions to the cleavage plane as evidence that the double 
chain ladder model is magnetically isotropic along the chain direction. In 
the strong coupling case two components of the triplet never contribute to 
the ground state at zero temperature, due to the strong single component 
contribution along the rungs. It has been suggested^^ that the triplet exci- 
tation can be considered as an analogue of Bose-Einstein condensation for 
magnons^*"^-'^ for this class of compoimds. The strongly coupled spin lad- 
der with magnon excitations for strong magnetic fields can be mapped to 
a one-dimensional XXZ-Heisenberg chain with an effective magnetic field. 
In this case the TBA equations reduce to only one level. The experimental 
magnetization curves'*'^ suggest an energy gap A « 31. IK and the critical 
field values Hd w 20 T and Hd » 50 T at T = 1.3 K. Fitting the zero tem- 
perature TBA critical fields and susceptibility to the experimental curves^ 
gives the coupling constants Jy = 5.5 K and J± = 57K for the integrable 
spin ladder model (1). 




(4) 



For the magnetic field along the x-direction they are given by 




(5) 
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Figure 2. Comparison between theoretical and experimental susceptibility curves versus 
temperature for the compound KCuCls. Circles and diamonds denote the experimental 
data extracted from Ref. 1 for an external field perpendicular or parallel to the chain 
direction. The solid and dashed curves are the corresponding susceptibility curves eval- 
uated directly from the HTE at H = OT. Fitting results in the coupling constants 
Jx = 57K and = 5.5 K with g = 2.29 (perpendicular), g = 2.05 (parallel) and 
fiB = 0.672 K/T. The conversion constant is Xhte ~ 0.40615xexp fixed in Ref. 15. 



3.1. Susceptibility 

The application of the HTE (3) for the free energy of Hamihonian (1) in- 
dicates that the couphng constants J\\ = 5.5 K and = 57K also give 
excellent fits to the susceptibility. The temperature dependence of the sus- 
ceptibihty curves is shown in Fig. 2. The solid and dashed lines denote 
the susceptibility for the external field perpendicular and parallel to the 
double chain direction, as derived from the free energy expression (3) with 
up to fifth order HTE. Here the Lande factors g — 2.29 (perpendicular) and 
g — 2.05 (parallel) for the external field direction were used. A rounded 
peak at T — 28.5 K in the zero magnetic field susceptibility curve indi- 
cates typical antiferromagnetic behaviour. The overall agreement with the 
experimental susceptibility curves is excellent. The susceptibility for the 
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magnetic field H(T) 

Figure 3. Magnetization versus magnetic field for the compound KCuCla with the same 
constants as in Fig. 2. This curve indicates the nature of the high field quantum phase 
diagram. The stiffness in the vicinities of the critical fields Hd and Hc2 is softened by 
increasing temperature. The critical fields predicted by the TBA are in good agreement 
with the experimental values. 



external field parallel to the chain direction has been examined via differ- 
ent theoretical models.^ Their conclusion favours a dimerized Heisenberg 
ladder structure with additional diagonal spin interactions, with the sug- 
gested coupling constants Jy = Jd = 8.35 K and J± = 50.1 K for the double 
chain structure compound. However, their fitting constants result in an en- 
ergy gap A w 38 K, which is much larger than the experimental value. 
We conclude that it is not necessary to introduce diagonal spin exchange 
interaction due to the strong dimerization along the rungs. The diagonal 
spin exchange interaction has only a weak effect on the low temperature 
behaviour. Moreover, the leg interaction is also suppressed by the relatively 
strong rung dimerization. 
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Figure 4. Specific heat versus temperature for different magnetic field strengths for the 
compound KCuCla with the same set of coupUng constants Jj^, Jy and g = 2.29. The 
sohd and dashed curves are evaluated from the HTE for H = OT and H = lOT. The 
conversion constant is Chte ~ 4.515Cexp- 



3.2. Magnetization 

The magnetization is a particular interesting quantity to study as the field 
dependent magnetization curve leads to the prediction of the low temper- 
ature phase diagram as well as magnetization plateaux. The high field 
magnetization curve evaluated from the TBA at zero temperature is shown 
in Fig. 3. By the nature of the high temperature expansion, we are unable 
to produce these very low temperature, T < 5.5 K, magnetization curves 
from the free energy (3) for this particular compound. This highlights the 
complementary role of the TBA and HTE approaches. The magnetization 
curve indicates that the rung singlets form a nonmagnetic ground state if 
the magnetic field is less than the critical field value Hd = 22.74 T. The gap 
closes at this critical point. If the magnetic field is above the critical point, 
the lower component of the triplet becomes involved in the ground state. 
The magnetization increases almost linearly with the field towards the criti- 
cal point Hc2 — 51.34 T, at which the ground state becomes fully polarized. 
This is in good agreement with the experimental values Hd « 20 T and 
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3.3. Specific Heat 

Fig. 4 shows the specific heat curves obtained from the HTE for the free 
energy at different magnetic field strengths. In the absence of a magnetic 
field the rounded peak indicates short range ordering with a large gap. At 
temperatures less than T = 17K the exponential decay signals an ordered 
phase. The magnetic field is seen to only weakly affect the magnetic spe- 
cific heat at low temperatures, mainly because of the strength of the rung 
singlets. As yet there appears to be no experimental data for the specific 
heat. 

4. Conclusions 

We have examined the magnetization, susceptibility and critical fields of 
the double chain compound KCuCls via the integrable spin ladder model 
(1). The theoretical results obtained from Thermodynamic Bethe Ansatz 
and High Temperature Expansion calculations are seen to lead to good 
agreement with the experimental measurements for these quantities. We 
conclude that this compound exhibits strong rung coupling which loads to 
dimerized rung spins. This is consistent with the experimental analysis. ^'^ 
We have also presented the specific heat curves for different magnetic fields. 
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